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ABSTRACT 

PSR B1259-63/SS 2883 is a binary system in which a 48-ms pulsar orbits around a 
Be star in a high eccentric orbit with a long orbital period of about 3.4 yr. Extensive 
broadband observational data are available for this system from radio band to very 
high energy (VHE) range. The multi-frequency emission is unpulsed and nonthermal, 
and is generally thought to be related to the relativistic electrons accelerated from 
the interaction between the pulsar wind and the stellar wind, where X-ray emission 
is from the synchrotron process and the VHE emission is from the inverse Compton 
(IC) scattering process. Here a shocked wind model with variation of the magneti- 
zation parameter a is developed for explaining the observations. By choosing proper 
parameters, our model could reproduce two-peak profile in X-ray and TeV light curves. 
The effect of the disk exhibits an emission and an absorption components in the X-ray 
and TeV bands respectively. We suggest that some GeV flares will be produced by 
Doppler boosting the synchrotron spectrum. This model can possibly be used and be 
checked in other similar systems such as LS I+61°303 and LS 5039. 

Key words: pulsars: individual: PSR B1259-63 — X-rays: binaries 



1 INTRODUCTION 

The discovery of the PSR B1259-63/SS 2883 system was first 
reported in 1992 (Johnston et al. 1992), and it is a binary 
system containing a rapidly rotating pulsar PSR B1259-63 
in orbit around a massive Be star companion SS 2883. The 
spin period of the pulsar is P = 47.76 ms and the spin- 
down luminosity is Lsp — 8 X 10^^ ergs s"\ The distance 
between the system and Earth is estimated to be about 1.5 
kpc (Johnston et al. 1994). 

This system is special for having unpulsed and non- 
thermal multi-band observational data spanning from radio 
band to TeV range, which are variable along with the or- 
bital phase (Johnston et al. 1994, 1996, 2005; Chernyakova 
et al. 2006, 2009; Uchiyama et al. 2009; Aharonian et al. 
2005, 2009). The broadband emission from this system is 
generally believed to be produced by the ultrarelativistic 
electrons accelerated in the shock from the interaction be- 
tween the pulsar wind and the outflow of Be star (Tavani 
& Arons 1997; Kirk, Ball & Skjasraasen 1999; Khangulyan 
et al. 2007; Takata & Taam 2009; Kerschhaggl 2010). Mas- 
sive stars usually produce very strong stellar wind during 
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their life time. For a Be star, the stellar wind is composed 
of two parts: a polar component and an equatorial compo- 
nent. The polar wind is faster and thinner; the equatorial 
component is slower and denser, and forms a disk surround- 
ing the Be star. The interaction between the pulsar wind 
and the stellar outflow will form a termination shock at the 
position where the dynamical pressures of the pulsar wind 
and the stellar wind are in balance, and this shock will ac- 
celerate electrons to relativistic velocities. These accelerated 
electrons will emit nonthermal emission via synchrotron pro- 
cess for X-ray regime (Tavani & Arons 1997) or IC [either 
synchrotron self-Compton (SSC) or external inverse Comp- 
ton scattering the thermal photons from the Be star (EIC)] 
process for VHE regime (Kirk, Ball & Skjaeraasen 1999). 
So this system is an important astrophysical laboratory for 
studying the interaction between the pulsar wind and the 
stellar wind, and the nonthermal radiation arising from the 
shocked pulsar wind. 

However, the prediction from a simple wind interaction 
model is not consistent with the observations. The simple 
wind interaction model predicts that the X-ray light curve 
reaches a maximum in flux at periastron (Tavani & Arons 
1997), but the observed X-ray light curve has a two-peak 
profile and has a drop in photon flux around the periastron 
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(Kaspi et al. 1995; Hirayama et al. 1999; Chernyakova et 
al. 2006, 2009). The behavior of TeV light curve is similar 
to that in X-ray range (Aharonian et al. 2005, 2009). Some 
authors suggest that hadronic scenarios may take place in 
the system (Kawachi et al. 2004; Neronov & Chernyakova 
2007). In the hadronic model. The Uhra-high energy (UHE) 
protons will be produced by the interaction between the pul- 
sar wind and the dense equatorial wind. These UHE protons 
can loose their energy in interactions with the protons from 
the stellar wind and produce pions. The decay of neutral 
pions will produce the TeV emission, and the leptons from 
the decay of charged pions will produce the radio to GeV 
emission via the synchrotron and IC processes respectively. 
In this situation, the two-peak profile in the light curve is 
related to the passage of the Be star disk. On the other 
hand, some other authors use some revised leptonic models 
to explain the drop of photon fiux towards periastron (Kan- 
gulyan et al. 2007; Takata & Taam 2009). Kangulyan et al. 
(2007) suggest that the observed TeV light curve can be ex- 
plained (i) by increasing of the adiabatic loss rate close to 
periastron or (ii) by the early sub- TeV cut-offs in the energy 
spectra of electrons due to the enhanced rate of Compton 
losses near the periastron. More recently, Takata & Taam 
(2009) explain the observed light curves by varying the mi- 
crophysical parameters in two different ways: (i) they use 
a given power-law spectral index p of electron distribution, 
and vary the magnetization parameter a and the pulsar wind 
bulk Lorentz factor F; (ii) the magnetization parameter and 
the power law index are varied for a given bulk Lorentz fac- 
tor. 

Among all the mechanisms, we believe that the varia- 
tion of microphysics model needs to be paid particular atten- 
tion, especially the variation of the magnetization parameter 
a. The varying of the magnetization parameter a along with 
the distance from the pulsar is a very natural hypothesis. In 
the Crab Nebula, a ~ 0.003 at a distance of ra ~ 3 x 10^^ cm 
(Kennel & Coroniti 1984a, 1984b). But at the light cylinder 
of a pulsar, a should be much larger and is about 10'' - 10^. 
Between the light cylinder and the termination shock, the 
magnetic energy will be gradually converted into the par- 
ticle kinetic energy, and a will vary along with the energy 
transition. Especially, a larger a around periastron will de- 
crease the energy in particles and make the cooling faster, 
which may cause a drop in photon flux. 

The PSR B1259-63/SS 2883 system is also special for 
having a large eccentricity of e = 0.87 and a long orbital pe- 
riod of about 3.4 yr. The semimajor axis of the orbit is about 
5 AU. Although the wind interaction model for this binary 
system is very similar to the pulsar wind nebulae (PWN) 
model (Kennel & Coroniti 1984a, 1984b), the detailed sit- 
uation is different. For the PWN around an isolated pulsar 
like Crab, the distance between the pulsar and the termina- 
tion shock Ta is about 0.1 pc and the magnetic parameter 
a is about 0.003 (Kennel & Coroniti 1984a, 1984b). How- 
ever, in the case of PSR B1259-63/SS 2883 system, due to 
the highly eccentric orbit, is in the range of 0.1 - 1 AU 
and cr should be accordingly larger. This gives us the unique 
opportunity to study the physical properties of the pulsar 
wind near the light cylinder where a is much higher. 

In this paper we develop a simple shocked wind model, 
in which the magnetization parameter a varies with rs to 
reproduce the broadband observational data. Besides PSR 



B1259-63/SS 2883, there are also some other similar binary 
systems being found for emitting the X-ray and VHE emis- 
sion, such as LS I-|-61°303 (Albert et al. 2006) and LS 5039 
(Aharonian et al. 2006). We suggest that this model may 
also be used in these similar systems. The outline of our 
paper is as follows: in Section 2, we introduce our model in 
detail. We then present our results and the comparison with 
observations in Section 3. Our discussion and conclusion are 
presented in Section 4. 



2 MODEL DESCRIPTION 

In our model, the broadband emission of PSR B1259-63/SS 
2883 system is mainly from the shock-accelerated electrons. 
Due to the interaction between the pulsar wind and the stel- 
lar wind, strong shocks will be formed, and the electrons will 
be accelerated at the shock front of the pulsar wind. This 
shock will also compress the magnetic field in the pulsar 
wind. These shocked relativistic electrons move in the mag- 
netic field and the photon fields (either from the synchrotron 
emission or from the Be star), and emit synchrotron and IC 
radiation to produce the multi-band emission. 

2.1 Stellar and pulsar winds, and termination 
shock 

The companion star of PSR B1259-63, which is named SS 
2883, is a massive Be type star. Be stars are characterized 
by strong mass outfiows around the equatorial plane where 
slow and dense disks will be formed (Waters et al. 1988). The 
velocity profile of the equatorial component can be described 
as 

R 

Ww,disk(-R) = U0,disk(-J5-)'", (1) 

it* 

where R is the distance from the stellar surface, «o,disk is the 
wind velocity at the stellar surface, -R* is the Be star radius, 
m is the outflow exponent. For a typical equatorial outflow, 
Wo, disk ~ 10^ cm s~', m is in the range < m < 2. On the 
other hand, there is also a polar wind component in the Be 
star outflow. Compared with the equatorial component, the 
polar wind has much higher velocity but lower density. The 
velocity distribution of the polar wind can be approximated 
as (Waters et al. 1988) 

Ww,polar(-R) = fO,polsar + {t-'oa — fO,polar)(l — R, / R)^ , (2) 

where Vaa is the terminal velocity of the wind at inflnity 
and VoD ~ 100«o,poisar, the index /3 is a free parameter. For 
a typical Be star, Voa ~ 10* cm s~^, and the value of /? is 
taken to be 1.5. 

The dynamical pressure of the stellar wind can be de- 
scribed as 

PAR) = Pv.{R)v„{Rf, (3) 

where pv,{R) = M /A-Kf„R?Vn]{R) is the wind density, 
M is the mass-loss rate of the Be star, /w is the out- 
flow fraction in units of A-k sr. For a typical polar wind, 
Mpoiar//w, polar ~ 10"** Mq yr"'; aud for a typical equa- 
torial wind, Mdisk//w,disk ~ 10"^ Mq yr"'. We assumed 
in our model that the density and velocity of the equato- 
rial wind evolve with the height above the equatorial plane 
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in Gaussian forms to match the polar wind. The dynamical 
pressure of the pulsar wind can be described as 



^'pui(r) = 



(4) 



where r is the distance from the pulsar, Lsp is the spin-down 
luminosity of the pulsar, /p is the pulsar wind fraction in 
units of 47r sr and c is the speed of light. We assume that 
the pulsar wind is isotropic and /p = 1 in our calculations. 
The location of the shock is determined by the dynamical 
pressure balance between the pulsar wind and the stellar 
wind, i.e.. 



= /9w(-Rs)|Uw(iis) - "orbl^, 

Alters 



(5) 



where is the distance between the termination shock and 
the pulsar, = d — is the distance of the shock from the 
stellar surface, d is the separation between the pulsar and 
its companion, Vorh is the pulsar orbital velocity (Tavani & 
Arons 1997). The value and the direction of Uorb are changed 
along with the orbital phase. 



2.2 Magnetic field 

In the preshocked pulsar wind at r^, the magnetic field Bi 
and the proper electron number density ni can be described 
as (Kennel & Coroniti 1984a, 1984b) 



a evolves with rg. The profile of the variation used in our 
model can be described as a power-law form. 



o- = crL( — ) 



(9) 



If we use the parameters of the Crab pulsar, i.e. Bl ~ 10® 
G and iVci ~ 10^'* s"\ which give ctl ~ 10^ at = 10*^ 
cm, and a — 0.003 at rs = 3 x lO^'^ cm, then the typical 
value of Qct is of order of unity. Note that here we obtain 
Oct by comparing with the Crab pulsar and assuming the 
conversion efficiency from the magnetic energy to the par- 
ticle kinetic energy is a simple power-law. But the actual 
condition in the PSR B1259-63/SS 2883 system may not be 
the same and larger or smaller values of may also be 
possible. So the exact value of is defined as a parameter 
in our calculations. According to the energy conservation, 
(1 + (j)Ti ~ constant, we define To as the value of Fi when 
(7 = and take it as a free parameter. We can see that Fi 
also evolves with rs due to the energy conservation. 

The downstream magnetic field is described as (Kennel 
& Coroniti 1984a, 1984b) 



B2 = Si W 1 



2 

U2 = 



1 



(10) 



80-2 + 1Q(J + 1 [64a2(a + 1)^ + 20o-(a + 1) + 1] 2 



16(cr + 1) 



(11) 



Bl 



ni = 



r|c(l -I- a) ' 



47ruiFir|mcC^(l -I- a) 



where Ui and Fi 



(6) 
(7) 

are the dimensionless ra- 
dial four velocity and the bulk Lorentz factor of the un- 
shocked pulsar wind, ttic is the rest mass of the electron, 
a = B\/{^'KT'\Uin\mcC^) is the magnetization parameter 
and is defined by the ratio of magnetic energy density and 
particle kinetic energy density in the pulsar wind. 

In the situation of the Crab Nebula, Kennel & Coroniti 
(1984a, 1984b) took the magnetization parameter a ~ 0.003 
at a distance of rs ~ 3 x 10^^ cm from the pulsar. In Tavani 
& Arons (1997), the magnetization parameter was adopted 
as cr = 0.02 for the PSR B1259-63/SS 2883 system. But at 
the light cylinder, the typical value of ctl should be as large 
as 



CTL = 



1.38X10^(^)^( 



2.3 Electron distribution 

It is usually assumed that the unshocked cold electron pairs 
can be accelerated to a power-law distribution in the ter- 
mination shock front, and be injected into the downstream 
post-shock fiow. So we adopt a power-law electron injection 
spectrum Q(7c) ~ (7c - 1)^'' (7c,min < 7c < 7e,max) in our 
work, where p is the electron distribution index and usu- 
ally varies between 2 and 3. Since the shock acceleration is 
a highly nonlinear process, it is very difficult to know how 
p evolves in different orbital phase. For simplicity, we as- 
sume that p is a constant. The minimum Lorentz factor can 
be determined from the conservations of the total electron 



number Lsp/Fme 



108cm 404'^ ^5.26 X lO^is-i'^ ' 



where Bl is the magnetic field at the light cylinder, rL 
is the radius of the light cylinder, N^± = NmNoj, Nm 



is the multiplicity and A^g.i 



5.26 X 10'^^(B/3 X 



10"G)(P/47.762ms)"^s"^ is the Goldreich- Julian particle 
flow at the light cylinder. In outer gap models (e.g. Cheng, 
Ho & Ruderman 1986a, 1986b; Zhang & Cheng 1997; 
Takata, Wang & Cheng 2010), the multiplicity due to var- 
ious pair-creation processes could reach 10** — 10^. We can 
imagine that between the light cylinder and the termination 
shock, the magnetic energy will be gradually converted into 
the particle kinetic energy. So in our work, we assume that 



/ Q(7c)d7c and the total electrons 



energy Lsp = / Q(7c)7c?ncC d7c (Kirk, Ball & Skjaraasen 
1999), and we can acquire 



7c, mill 



Fi 



p~l 



(12) 



The maximum Lorentz factor of electrons can be determined 
by equating the cooling timescale of electrons fcooi with the 
particle acceleration timescale tac as the following form. 



7c, max 



67re 



cttB2(1 + K)' 



(13) 



where e is the electron charge, ctt is the Thompson scat- 
tering cross section, Y is the Compton parameter and is 
defined below. The cooling timescale fcooi ~ 7c'TicC^/-Prad, 
where P^ad = (4/3)CTTC7„^(B|/87r)(l -|- Y) is the radia- 
tion power of electrons. The particle acceleration timescale 

tac — 7e'TleC/ei32. 

The cooling of electrons due to synchrotron and IC ra- 
diation will steepen the distribution of electrons above a 
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critical Lorentz factor 70, c, which can be derived from equat- 
ing the cooling timescale tcooi with the dynamic flow time 
Tdyn of electrons. So 70, c can be expressed as (Sari, Piran & 
Narayan 1998) 

(TT-B2^dyn(l + Y) 

In our model, we assume the dynamic flow time of electrons 
in the radiation cavity is Tdyn = ^rg/wt, where ^ is a di- 
mensionless parameter and Uf is the postshock flow velocity. 
In Tavani & Arons (1997), they used ^ = 3 by taking into 
account the nonspherical shape of the shocked region and 
Vi = c/3. The electrons with Lorentz factor 7c > 70, c will be 
cooled by the radiation rapidly. 

The Compton parameter Y used above describes the 
effect of IC (Sari & Esin 2001), which is deflned as the ra- 
tio of the IC luminosity (including SSC and EIC) to the 
synchrotron luminosity, 

■^r, N isSC + -f/EIC 
y(7e) = J 

_ Usyn[v ^ I^Kn(7c)] + UstB.iW < !^Kn(7c)] 



Usynif ^ ^T<n(7o)] = 
Usts.vV' < l'KN(7e)] = 



UB 

1 + y(7e) iTTCrl ' 
Lsta^ilu ^ ^'Kn(7c)] 

47rc(d — rs)2 



Ub 



D2_ 

Stt' 



(15) 
(16) 
(17) 
(18) 



where i/KN(7c) ~ mcC^/{h'yc) is the critical frequency of 
scattering photons above which the scattering with electrons 
of Lorentz factor 7e enter the Klein-Nishina (KN) regime, 
vile) = 7?rad??KN(7o) IS the radiatlou efficiency, 7?rad is the 
fraction of the electron's energy which is radiated away, 
and 77kn(7o) is the fraction of synchrotron photons below 
the KN limit frequency i'KN(7e). In the slow cooling case, 

7e,min < 7c,c, ??rad = (7o,c/7e,miii) (WaUg et al. 2010), 



In the fast cooling case, 7e 
2010), 



)3 , !^Kn(7o) < fm, 

< !^Kn(7c); 
< 7e,min ,r?rad = 1 (WaUg Ct al. 



i^t^i^^t'^ /^kn(7o)^^c, 

'7KN(7c)=<i C-^T^t^. /^c < //KN(7e) < i^,n, (20) 

1, Vvn ^ !^Kn(7c), 

where and are the typical synchrotron radiation fre- 
quencies of electrons whose Lorentz factors are 7c,min and 
7e,c respectively, and are deflned as Vm ~ 37c,me-B/(47rmcc) 
and Vc = 37e ce5/(47rmec) respectively. 

The radiation revised electron spectrum N{'yc) can be 
obtained from the continuity equation of the electron distri- 
bution (Ginzburg & Syrovatshii 1964), 



dn{y^,t) _^ a7en(7c,t) ^ ^^^^^ 



(21) 



dt djc 

Similar to the LS I-|-61°303 system, because the cooling 
and dynamic flow timescales are much smaller than the or- 
bital period in the PSR B1259-63/SS 2883 system, we use 



dn{yc,t)/dt = to calculate the electron distribution at a 
steady state (Zabalza, Parades & Bosch-Ramon 2011) and 
obtain, 

(1) 7c, c ^ 7c,min ^ 7G,max, 



iV(7c) oc 



(7c - 1)- 



7c, c < 7c < 7o,min, 

"'e,max ; 



(7c - 1) '■^'^^^ , 7c,min < 7o < 7' 



(22) 



]\r/ \ „, ] (7c 1) ^ I 7c,min ^ 7c ^ 7o,c, ,„q\ 



(3) 'yc,c ^ 'yc.,m.a.x.-) 

iV(7c) (X (7c - 1)-^ 7 '7c,min «; 7c 5? 7c 



(24) 



The coefficient can be calculated from A^'tot = / Af(7c)d7c, 
where iVtot = (Ari/47r)Lsprdyn/[rimcC^(l-|-(j)], and Ari/47r 
is the fraction of pulsar wind electrons which are acceler- 
ated to relativistic velocities in the shock front. In our cal- 
culations, we assume that the typical scale of the radiation 
cavity is Vs for simplicity, and therefore AQ ~ 7r/4. 

2.4 Radiation process 

In our model, the multi-band photons of PSR B1259-63/SS 
2883 are from the synchrotron radiation and IC radiation 
(including SSC scattering of the synchrotron photons and 
EIC scattering of the thermal photons from the Be star) of 
the shock-accelerated electrons. 

The synchrotron radiation power at frequency v from a 
single electron with Lorentz factor 7c is given by (Rybicki 
& Lightman 1979) 



,2 



(25) 



where Uc — 37cei3/(47rmcc). The function F{x) is defined 



F(x) = X 



A'5/3(fc)dfc. 



(26) 



The SSC radiation power at frequency v from a single elec- 
tron with Lorentz factor 7c is given by (Blumenthal & Gould 
1970) 



^'.^'^'^(7e) = 3aT 



di^s . „ g{x,y), 



g{x, y) = 2ylny + {1 + 2y){l - y) + 



2 2 

X y 



2(1 + xy) 



(27) 



(1-y), (28) 



where f^J^ = / ^^7'~'(7c)A^(7c)d7c/47rrJ is the flux den- 
sity of the synchrotron radiation, x — ^'^chva/ ijn^c?), y = 



hu/[x{'~icmcC - hu)], Va.min = fnieC /[47e(7e 



Correspondingly, the EIC radiation power at frequency u 
from a single electron with Lorentz factor 7c is given by 
(Aharonian & Atoyan 1981; He et al. 2009) 



J v. 

h{i,he) = 1 + 



.STAR 

2i 2^2 



(29) 
(30) 



2(1-0 6^(1-0^' 
where jl'^'^^' = TTB^{Tca)(R*/Rsf is the flux density of the 
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Be star photons, Bi,{Tcs) = 2hv^ / [exp{hi' / kTeg) — 1] is 
the brightness on the Be star surface, Tog is the effective 
temperature of the star, h is the Planck constant, k is the 
Boltzmann constant, ctt is the Thompson cross-section, ^ — 
hv / {'yemc(?) , be = 2(1 — cos9sc)'yehi'a /{meC^), hvs ^ hu ^ 
7emeC^&9/(l + bg), Bsc is the angle between the injecting 
photons and the scattered photons, and is varied along with 
the orbital phase. Note that here the Be star is assumed as 
a point-like and black body emitter for simplicity. 

The observed total flux densities from all the shock ac- 
celerated electrons at frequency v is 



_F'^(cos6'sc) = 



/'7e,ma 
J min{'ye 



ni„) 



d7eiV(7e)[Pr''(7a 



+ f'f '^(le) + '^(7e,COS0Sc)]. 



(31) 



As the pulsar orbiting around the Be star, the spatial condi- 
tion and the parameter cos^sc will be changed. Some phys- 
ical conditions, such as the the electron distribution N{-ya), 
the downstream magnetic field and the photon field, will 
vary accordingly. So a variable flux along with the orbital 
phase will appear. 

The gamma-ray photons produced by the IC process 
could be absorbed in the dense stellar photon field through 
pair production (Gould & Schreder 1967). The opacity could 



be roughly estimated by 



a-yjUtRis ~ 0.67 at perias- 



tron, where a^-y « tt/5, n« is the stellar photon density. 
Because the opacity at the periastron is already the largest 
value in all orbital phases but it is still less than unity, we 
will ignore the effect of pair production in our calculations 
for simplicity. 
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3 RESULTS 

In this section, we will present some calculated results using 
our model, and compare them with the observations. Some 
parameters used in our calculation for the PSR B1259-63/SS 
2883 system are as follows (Tavani & Arons 1997); For the 
orbital parameters, we take the eccentricity e = 0.87, the 
semimajor axis a = 4.8 AU; For the compact object PSR 
B1259-63, we take the spin-down luminosity Lap = 8 x 10^^ 
erg s^^; For the Be star SS 2883, we take the stellar lumi- 
nosity Latar = 5.8 X 10*^0, the stellar mass Matar = IOM0, 
the stellar radius _R* = IORq, the effective temperature of 
the star Tad = 27000 K. The distance between the system 
and Earth is taken to be 1.5 kpc (Johnston et al. 1994). We 
use ^ = 3 and v{ — c/3 in our calculations, which are the 
same as those suggested by Tavani & Arons (1997). 

3.1 Variation of some shock parameters 

We will first give some results on the variations of the dis- 
tance between the shock and the pulsar Va, the magnetiza- 
tion parameter a and the electron number A'^i_iokoV whose 
typical synchrotron emission frequencies are in the range of 
1-10 keV as a function of orbital phase, as shown in Figure 1. 
The periastron is taken to be at orbital phase ~ 0° through- 
out the paper. In solid lines, we only consider the polar wind 
and the input parameters we used are as follows: the termi- 
nal velocity of the polar wind — 1.3 x 10* cm s~^, the 
index of the polar wind /3 = 1.5, Afpoiar//w, polar = 2 x 10~* 



Figure 1. Top panel: The variation of shock distance along 
with the orbital phase; Middle panel: The variation of magneti- 
zation parameter a along with the orbital phase; Bottom panel: 
The variation of electron number A^i^iokcV along with the orbital 
phase. The solid and dotted lines correspond to the conditions 
without and with the effect of disk respectively. 

Mq yr~^, the electron distribution index p — 2.2, the 
magnetization parameter at light cylinder (rg = lO^cm) 
o-L = 8 X 10^, the index = 1.0 and To = 2 x 10^ In 
dotted lines, we add an equatorial component with the ini- 
tial velocity uo.disk ~ W^' cm s^^, the index m — 1.5 and 
Mdiak//w,diBk = 2 X 10-^ Mq yr-\ In our calculations, the 
half-opening angle of the disk (projected on the pulsar or- 
bital plane) is A6'diak = 10°, and the intersection between 
the stellar equatorial plane and the orbital plane is inclined 
at 6'disk ~ 75° to the major axis of the pulsar orbit, which 
are in the range suggested by Chernyakova et al. (2006) with 
A^diak — 18°. 5 and 6'di3k — 70°. These input parameters we 
used here are chosen by modeling the observations (See be- 
low) . It can be seen that because the shock distance varies 
along with the orbital phase, the magnetization parameter 
(7 and the electron number A'^i_iokoV evolve accordingly. We 
can also see from the bottom panel of Figure 1 that because 
the total electron number evolves as A'tot oc Vs, A'^i-iokcV 
varies as Ai'i_iokcV oc when ra is large. As the pulsar ap- 
proaching the periastron, the IC process is more and more 
important and the value of y(7e,c) will increase. At a cer- 
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Figure 2. The calculated spectra. The solid line and the dashed 
line correspond to the conditions at periastron and at the orbital 
phase of 100° respectively, where the effect of disk is not consid- 
ered; the dotted line corresponds to the result in the disk. 

tain position, some electrons will enter the fast cooling case 
(7e,c < 7e,min) with the distribution index being — (p + 1) 
instead of —p in the slow cooling case. At this moment, the 
electrons will be distributed in a more broad range, and the 
value of A'^i-iokoV will decrease more rapidly around the pe- 
riastron because more and more electrons are cooled to the 
lower energies. In the dense disk, and a reach a mini- 
mum and a maximum respectively, and A'^i-iokcV reaches a 
minimum because of the smallest electron number and the 
fastest cooling effect. 

3.2 Spectrum 

The calculated spectra from 10~^ eV to 100 TeV are illus- 
trated in Figure 2. The input parameters used here are the 
same as those in Figure 1. The different types of lines repre- 
sent the spectra in different orbital phases: the solid line and 
the dashed line correspond to the conditions at periastron 
and at the orbital phase of 100° respectively, where the effect 
of the disk is not considered; the dotted line corresponds to 
the result in the disk (the orbital phase = —75°). 

The emission below 100 MeV is mainly from the syn- 
chrotron radiation. The shape of synchrotron spectrum is 
mainly determined by the magnetic field and the electron 
distribution. It can be seen that between 1-10 keV the flux 
in the dotted line is the highest in all the three lines. Al- 
though the electron number A'^i-iokeV reaches a minimum 
in the disk, the radiation flux can be very large because the 
magnetic field is highest at this moment, as shown in Figure 
1. The large flux in the disk will produce an "emission" fa 
profile in the X-ray hght curve. We can also see that in solid 
line the spectrum below 10 eV is softer and the flux below 1 
eV is higher than those in other two lines, which is because 
more electrons are cooled to this energy range at periastron. 

The emission above 1 GeV is mainly produced by the 
EIC process. The calculated SSC flux is lower than the EIC 
flux by more than three orders of magnitude for the rea- 
son that the SSC emission is strongly suppressed by the 
KN effect, so we neglect the effect of SSC in our calcula- 
tions throughout the paper. There are two breaks in the EIC 



spectra corresponding to the two critical energies 7e,min?TieC^ 
and {rrieC^Y /kTefi (Yu 2009). In our calculations, 7e 
corresponds to the first break and (meC^)'^ /kTaS ~ 0.1 TeV 
corresponds to the second break in the EIC spectra. We 
can see that the flux in the dotted line is the lowest in all 
the three lines. The reasons are as follows: (1) In the disk, 
the magnetization parameter a gets a maximum and 7o,min 
reaches a minimum accordingly. The small 7G,min makes the 
EIC spectrum begin to break at a low energy; (2) the elec- 
tron number whose EIC emission frequencies are above 1 
GeV decreases duo to the rapid cooling in the disk; (3) The 
strong wind in the disk pushes the shock to the pulsar, so the 
photon density from Be star decrease. The low flux in the 
disk will produce an "absorption" profile in the TeV hght 
curve. 

The Fermi Gamma-ray Space Telescope observed the 
PSR B1259-63/SS 2883 system recently (Tam et al. 2010; 
Abdo et al. 2010a, 2010b, 2011b; Kong et al. 2011). The 
Large Area Telescope (LAT) onboard Fermi covers the en- 
ergy band from 20 MeV to greater than 300 GeV. We can 
see that in our model, the emission in this energy range is 
produced by the synchrotron and the EIC processes together 
(See Section 3.6 for detailed discussion). 

3.3 Comparison with observations in the 1-10 keV 
band 

The X-rays from the PSR B1259-63/SS 2883 system was first 
detected by ROSAT (Cominsky, Roberts & Johnston 1994). 
After that, the ASCA satellite also observed this system in 
X-rays (Kaspi et al. 1995; Hirayama et al. 1999). The re- 
sults from the ASCA satellite show that the X-ray flux from 
the source is highly variable in different orbital phase, and 
the hght curve has two peaks before and after the perias- 
tron respectively. More recently, new X-ray observation data 
of this source were reported by Chernyakova et al. (2006, 
2009). The results include the XMM— Newton observations 
at the beginning of 2004, some earlier unpublished X-ray 
data from BeppoSAX and XMM — Newton (Chernyakova 
et al. 2006), and the unprecedented detailed observations 
by Suzaku, Swift, XMM — Newton and Chandra mis- 
sions during the 2007 periastron passage (Chernyakova et 
al. 2009). These new data are well in agreement with pre- 
vious observations, and show that the orbital light curve 
in X-rays does not exhibit strong orbit-to-orbit variations 
(Chernyakova et al. 2009). 

Here, we use our model described in Section 2 to repro- 
duce the X-ray light curve of the PSR B1259-63/SS 2883 
system. We assume that the X-ray emission is mainly from 
the synchrotron radiation of the relativistic electrons. The 
X-ray data are taken from Chernyakova et al. (2006, 2009). 
In this subsection, we will only use the polar component as 
the stellar wind. The role of the disk will be discussed later. 
The calculated hght curves and the observed data in 1-10 
keV band are presented in Figure 3. Except for special dec- 
laration, the parameters we used here are the same as those 
in the solid lines of Figure 1. In the dotted line, we use a 
constant value of o" = 0.02 for reference purpose. Differently, 
in the solid and dashed lines, we make the magnetization pa- 
rameters vary as cr = o-i^{r /ri^)~°''' . We can see that when 
we use a constant magnetization factor a — 0.02, as sug- 
gested by Tavani & Arons (1997), the resulting hght curve 
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Figure 3. The calculated X-ray light curves and their comparison 
with observations. The solid line and the dashed line correspond 
to Oct = 1.0 and = 2.0 respectively. The dotted line corre- 
sponds to a constant magnetization parameter cr = 0.02. 

cannot explain the two major characteristics in observations: 
(1) The observed 1-10 keV flux has a sharp rise of about 20 
times from apastron to pre-periastron, but the ratio of flux 
between the maximum and the minimum is only about 5 
in the dotted Une; (2) In the dotted line, the flux get the 
maximum at periastron, in opposition to the observations 
where the flux has a minimum around the periastron. When 
the variation of cr is considered, the calculated results flt the 
above two characteristics much better. The appearance of 
this type of bimodal structure is due to the joint action of 
the variable magnetic held and electron number. As shown in 
Figure 1, as the pulsar approaching the periastron from the 
apastron, the magnetization parameter o becomes larger, so 
the magnetic field B and the synchrotron radiation flux in 
X-rays (-Lx) should increase more rapidly than those in the 
situation where a is a constant. A much sharper rise of about 
20 times in flux could be produced at this moment. On the 
other hand, a larger magnetic field B and a smaller shock 
distance around the periastron will make the cooling of 
particles stronger. Much more electrons will be cooled to the 
lower energies, and the electron number A''i_iokoV decreases 
more rapidly accordingly, as shown in Figure 1. So the flux 
in the 1-10 keV range has a minimum around the periastron. 

The solid and dashed lines in Figure 3 illustrate the 
effect of the index on the 1-10 keV hght curves. The 
parameters used in solid line are the same as those used 
in the solid line of Figure 1. We choose the value of the 
magnetization parameter tj as a constant at periastron, and 
take ao- = 2.0 in the dashed line. It is clearly seen that when 
Oct becomes larger, the rise and the drop of flux around the 
apastron become sharper. This is due to the more rapid 
change in the value of cr. 

3.4 Comparison with VHE observations above 1 
TeV band 

As early as 1999, Kirk, Ball & Skjaeraasen (1999) predicted 
that the PSR B1259-63/SS 2883 system could produce VHE 
radiation, and the light curve should reach the maximum 
around the periastron. About 5 years later, the TeV 7-rays 
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Figure 4. The calculated TeV light curves and their comparison 
with observations. The solid line and the dashed line correspond 
to Oct = 1.0 and = 2.0 respectively. 

from this source were first detected by the High Energy 
Stereoscopic System (HESS) around its periastron passage 
in 2004 (Aharonian et al. 2005). This makes the PSR B1259- 
63/SS 2883 system be the first known binary system to emit 
photons above 1 TeV (Aharonian et al. 2005). The observed 
TeV light curve is similar to that in 1-10 keV for having two 
peaks in pre- and post-periastron respectively with a mini- 
mum in fiux around the periastron (Aharonian et al. 2005), 
and is contrary to the prediction by Kirk, Ball & Skjaraasen 
(1999). The TeV light curve was also found to have signifi- 
cant variations on timescales of days, which makes the sys- 
tem be the first variable VHE source in our Galaxy (Aharo- 
nian et al. 2005). The HESS team also observed the source 
in 2005, 2006 and around the 2007 periastron passage (Aha- 
ronian et al. 2009). They found that there is no significant 
VHE 7-ray flux near apastron in the 2005 and 2006 obser- 
vations, and conflrmed that this system is a variable TeV 
emitter. 

Here, we use our model introduced in Section 2 to repro- 
duce the TeV light curve of the PSR B1259-63/SS 2883 sys- 
tem. We consider that the photons above 1 TeV are mainly 
produced by the EIC process. Our modelling light curve are 
presented in Figure 4 and the parameters used here are the 
same as those used in Figure 3. The observational TeV data 
are taken from Aharonian et al. (2005, 2009). We can see 
that our model can explain the two-peak profile in TeV ob- 
servations qualitatively. 

3.5 Effect of the equatorial component in the 
stellar wind 

Be stars like SS 2883 are characterized by strong mass out- 
flows around the equatorial plane. They can produce disks 
which are slower and denser than the polar wind (Waters et 
al. 1988). The existence of a disk in the PSR B1259-63/SS 
2883 system has been conflrmed by the radio observations 
(Johnston et al. 1996, 2005). The disk of Be star SS 2883 in 
this binary system is believed to be tilted with respect to the 
orbital plane. The line of intersection between the disk plane 
and the orbital plane is oriented at about 90° with respect to 
the major axis of the binary orbit, but the inclination of the 



8 S. W. Kong, Y. W. Yu, Y. F. Huang and K. S. Cheng 



disk is not constrained (Wex et al. 1998; Wang, Johnston & 
Manchester 2004). Chernyakova et al. (2006) analyzed the 
X-ray light curve of the system detailedly and suggested 
that the disk has a well-defined Gaussian-profile. The half- 
opening angle of the disk (projected on the pulsar orbital 
plane) is A^disk — 18°. 5, and the intersection between the 
stellar equatorial plane and the orbital plane is inclined at 
^disk — 70° to the major axis of the pulsar orbit. 

The presence of disk not only affects the radio radia- 
tion from the system, but also plays an important role in 
producing the X-ray and VlfE 7-ray emission (Chernyakova 
et al. 2006). Ball et al. (1999) suggested that the acceler- 
ation process will be most efficient when the pulsar passes 
through the disk. Kawachi et al. (2004) suggested that in 
addition to the relativistic electrons, VHE protons could be 
produced by the interaction between the pulsar wind and 
the dense equatorial wind. These hypothesis predict that 
there will be an increase of flux in disk in both the X-ray 
and TeV bands and can explain the two-peak profiles in the 
observed light curves. Recently, Kerschhaggl (2010) found 
that there is a drop in the TeV light curve when the pulsar 
passes through the dense disk. They suggested that this is 
due to an increasing nonradiative cooling effect in the disk. 

Here we add the equatorial component of the stellar 
wind in our calculations, and assume the mass-loss rate and 
velocity of the equatorial wind have Gaussian profiles to 
match those of the polar wind. Our results in the X-ray and 
TeV bands are presented in Figure 5. The solid and dotted 
lines correspond to the conditions with and without equato- 
rial components respectively. The parameters used here are 
the same as those used in Figure 1. We can see in Figure 
5 that the X-ray flux increases in the passage of disk, but 
the flux in TeV range decreases significantly at this moment, 
which is consistent with the analysis by Kerschhaggl (2010). 
As discussed in Section 3.2, the appearance of an emission 
component in the X-ray light curve is mainly due to the com- 
petition between the decrease of the electron number and the 
increase of the magnetic field inside the disk, whereas the 
appearance of an absorption component in the TeV range is 
due to the low break energy, the decrease of electron num- 
ber resulting from rapid cooling and the decrease of photon 
density in the disk. 

3.6 Observations by LAT onboard Fermi 

The PSR B1259-63/SS 2883 system passed through its peri- 
astron in mid-December 2010 recently. The LAT onboard 
Fermi observed this source detailedly. An average flux 
(above 100 MeV) of 1.2±0.6 x 10"^ ph cm"^ s"^ and a pho- 
ton index of 2.2 ±0.2 in the period 2010-11-17 to 2010-12-19 
UTC was reported by the Fermi team (Abdo et al. 2010b). 
They also reported a flux (above 100 MeV) of 1.7±0.4x 10"** 
ph cm~^ s~^ with a photon index of 3.0 ± 0.3 in the pe- 
riod 2011-01-17 to 2011-01-19 UTC, and a flux (above 100 
MeV) of 1.8 ± 0.2 x 10"® ph cm"^ s"^ with a photon in- 
dex of 3.1 ±0.2 in the period 2011-01-14 to 2011-01-19 UTC 
(Abdo et al. 2011b). Some other analyses of LAT data were 
also done. A flare-like profile with a fiux (300 MeV - 100 
GeV) of around 4 x 10~* ph cm~^ s~^ and a photon index 
of about 1.7 in the period 2010-11-18 00:00:00 (UT) to 2010- 
11-21 00:04:42 (UT) (Tam et al. 2010), and a fiux (200 MeV 
- 100 GeV) of 4.8 ± 0.9 x 10"^ ph cm"^ s"^ with a photon 
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Figure 5. The effect of disk on the X-ray and TeV light curves. 
The solid and dotted lines correspond to the conditions with and 
without the effect of disk respectively. 
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Figure 6. The calculated GeV light curves and their comparison 
with observations. The emission in solid line is from synchrotron 
radiation and the emission in dashed line is from the EIC process. 



index of 3.6 ± 0.4 in the period 2011-01-14 to 2011-01-16 
UTC (Kong et al. 2011) were reported. 

As shown in Figure 2, the emission between 100 MeV 
and 300 GeV range is produced by the synchrotron and the 
EIC processes together. In Figure 6, we present our calcu- 
lated light curve between 100 MeV and 300 GeV using the 
parameters which are the same as those in the dotted lines of 
Figure 1. The emission in solid line is from synchrotron ra- 
diation and the emission in dashed line is from the EIC pro- 
cess. The observation data are taken from The Astronomer's 
Telegrams (ATel; Tam et al. 2010; Abdo et al. 2010a, 2010b, 
2011b; Kong et al. 2011), and the fiux derived from Tam et 
al.'s (2010) and Kong et al.'s (2011) reports have been ex- 
trapolated to the 100 MeV and 300 GeV energy range using 
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photon indices of 1.7 and 3.6 respectively. We can see tiiat 
the calculated GeV flux are dominated by the EIC process, 
but the EIC flux is still lower than the observations. This 
may be due to two reasons: (1) In the ATel reports, a single 
power-law has been assumed for the GeV photon spectrum, 
and we use it to estimate the energy flux in the flgure. But 
in our calculation the spectrum is not a single power-law. 
It may overestimate the GeV flux by using a single power- 
law spectrum in the flgure. Some recent analyses actually 
show that the spectrum between 100 MeV and 300 GeV is 
not a single power-law (Adbo et al. 2011; Tarn et al. 2011). 
(2) Some observations with high flux may be flares, which 
cannot be produced in the steady model. In other words 
the flares should be emitting from different region of the 
shock, which is also suggested by Tam et al. (2011) com- 
paring the 100 MeV to 300 GeV spectra with the > 300 
GeV gamma-ray spectra. We suggest that the flares may 
come from the relativistic Doppler-boosting effect (Dubus, 
Cerutti & Henri 2010). Bogovalov et al. (2008) presented 
a hydrodynamic simulation of the interaction between rela- 
tivistic and nonrelativistic wind in PSR B1259-63/SS 2883 
system. They found that the electrons in the shocked pul- 
sar wind can be accelerated from bulk Lorentz factors ~ 1 
around the termination shock to nearly 100 far away. The 
spectrum of synchrotron emission in our calculations peaks 
at about 0.1 GeV. With such a large bulk Lorentz factor 
of ~ 100, the peak will be Doppler-boosted to about 10 
GeV, which could be detected by Fermi satellite. However, 
emission from the electrons with a bulk Lorentz factor of 
100 should be strongly beamed. When the line-of-sight is in 
the beaming direction, we can receive the GeV flux, other- 
wise, the GeV photons disappear. In this case, a flare pro- 
flle may be found. The much higher intensity of the flares 
is mainly due to the beaming effect because instead of ra- 
diating isotropically the GeV photons from the large bulk 
Lorentz factor region is highly beamed at a much smaller 
solid angle. If the geometry of the shock has a simple bow 
shock structure at the tail, we expect that two flares will 
occur when the line of sight passes through the two edges 
of the cone of the bow shock and each flare can last about 
an orbital phase 2/r ~ 12°(r/10)"\ where T is the bulk 
Lorentz factor at the tail of the shock. If the photon num- 
bers are sufficiently high during the flares we should be able 
to find the synchrotron cut-off energy (i?syn ~ 100 MeV cf. 
Figure 2) boosting to VEsyn ~ 1 GeV (F/IO) in these flares 
and the general spectrum in Fermi range should satisfy a 
simple power with exponential cut-off. 

Recently, Tam et al. (2011) found the flux in the flaring 
period in Fermi observations is enhanced by a factor of 5-10, 
which suggests a Doppler factor of around 1.5-2. With a low 
energy break around 10 keV observed by Suzaku (Uchiyama 
et al. 2009), the enhancement has a factor of only about 3 in 
the X-ray band, which is smaller than that in the GeV band 
and is less obvious. In addition, the flaring period is near 
the second peak in the X-ray light curve. This peak may 
be produced by the effect of Doppler boosting and disk to- 
gether. For the TeV flux, because more IC scattering occur 
in the KN regime for F > 1, the TeV flux should be sup- 
pressed. In addition, the stellar photon density in the bow 
shock tail is smaller than that in the shock apex, so the IC 
emission should be lower accordingly. The increase of TeV 
flux due to the Doppler boosting and the decrease due to 



the KN effect and smaller photon density will compensate 
each other, and the Doppler boosting may not modulate 
the TeV light curve signiflcantly. Furthermore, because the 
emission from the shock tail is temporally and spectrally 
different with that from the shock apex, the magnetic flelds 
and the electron distributions may not be the same in these 
two ranges. 



4 CONCLUSION AND DISCUSSION 

The PSR B1259-63/SS 2883 system is an attractive binary 
system, and is also an unique astrophysical laboratory for 
probing the time dependent interaction between the pulsar 
wind and stellar wind, and studying the physics of the pulsar 
wind in a different range compared with an isolated pulsar. 
In this paper, we have modeled the X-ray and TeV obser- 
vations of the source using the wind interaction model un- 
der the hypothesis that the X-ray emission is mainly from 
the synchrotron process and the TeV photons are mainly 
produced by the EIC effect where the relativistic electrons 
in the shock up-scatter the photons from the Be star. The 
effect of the disk exhibits an emission and an absorption 
components in the X-ray and TeV bands respectively. More 
importantly, we assume that the microphysics parameters 
are not constant and take the magnetic parameter a vary- 
ing as a power-law form along with the distance between 
the termination shock and the pulsar. This is a reasonable 
assumption because a has a very large value at the light 
cylinder and reaches a very small value of about 0.003 at a 
distance of about 0.1 pc from the pulsar (Kennel & Coro- 
niti 1984a, 1984b). The variation of microphysics parameters 
has also been suggested in some other astrophysical phenom- 
ena like Gamma-ray bursts (Kong et al. 2010). We also try 
to explain the GeV light curve which is observed by LAT 
very recently, and suggest that the actual photon spectra 
may not be power-law distributed as reported in the ATel 
and some flare-like proflle may come from the relativistic 
Doppler-boosting effect. 

In a recent paper, van Soelen & Meintjes (2011) inves- 
tigated the effect of the infrare emission from the Be star's 
disk on the inverse Compton process. The free-free and free- 
bound emission occurred in the disk is thought to produce an 
infrared excess in the Be star's photon spectrum. The scat- 
tering of this infrared excess will influence the gamma-ray 
production. As discussed by van Soelen & Meintjes (2011), 
for the electrons with Lorentz factors 70 > 10^, the scat- 
tering of photons around the peak in the stellar spectrum 
will be in the KN limit, while the scattering with photons 
in the infrared excess will be in the Thompson limit and in- 
crease the emission at GeV band. But for the broad electron 
distribution with 7o,min ~ 10^, the scattering with photons 
around the peak of the stellar spectrum will occur in the 
Thompson limit, and most of the GeV photons are produced 
by the Thompson scattering of the whole stellar spectrum. 
So the effect of the infrared excess on the spectrum could 
be negligible. In our calculations above, the electron distri- 
butions are very broad with min(7o,min, 7e,c) between a few 
X 10^ to a few x 10*. The resulting influence on the GeV flux 
by the infrared excess is insigniflcant, so we do not consider 
this effect in our calculations. 

There are also some other binary systems in our Galaxy 
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similar to PSR B1259-63/SS 2883 being found for emitting 
the X-ray and VHE emission, such as LS 5039 (Aharonian 
et al. 2006), LS I+61°303 (Albert et al. 2006) and possibly 
HESS J0632-I-057 (Hinton et al. 2009). Our model can possi- 
bly be used in explaining the broadband observational light 
curves in these systems. However, the orbital properties are 
different in these systems. In PSR B1259-63/SS 2883, the 
separation between the pulsar and the massive star is 0.7-10 
AU. But the separation is only about 0.1-0.2 AU in LS 5039 
and 0.1-0.7 AU in LS H-61°303 (Dubus 2006). The mag- 
netization parameter a should be larger accordingly. The 
detailed study of the radiation from the LS 5039 and the LS 
I-|-61°303 systems could help us to understand the physical 
condition of the pulsar wind in the position where is much 
closer to pulsar compared with the PSR B1259-63/SS 2883 
system. 
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